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ABSTRACT 

An  objective  technique  is  proposed  to  use  digitized  satellite  cloud 
brightness  data  to  estimate  large-scale  flow  patterns  over  data-void 
tropical  regions.   It  is  basically  a  diagnostic  model  which  assumes: 
1)  area-averaged  cloud  brightness  is  positively  correlated  with  large- 
scale  divergence  in  the  tropical  upper  troposphere;  and  2)  large-scale 
tropical  flow  is  quasi-barotropic  and  quasi-non-divergent.   It  can  be 
used  at  any  upper  level  where  divergence  is  important  in  determining 
vorticity.   The  model  requires:  1)  boundary  conditions  determined  from 
surrounding  wind  reports,  2)  a  mean  zonal  flow  determined  from  clima- 
tology, and  3)  an  equivalent  divergence  forcing  function  constructed 
empirically  from  the  brightness  data. 

The  technique  is  tested  daily  over  a  western  North  Pacific  region 
for  July-August  1971.   Results  show  that  for  25%  of  the  days  tested, 
the  model  produces  a  flow  field  which  accurately  resembles  the  major 
features  of  the  stream-function  field  analyzed  by  NMC.   In  another  30% 
of  the  days  it  provides  some  valuable  information  about  the  flow 
patterns  which  would  be  difficult  to  obtain  from  boundary  information 
alone.   Experiments  also  suggest  that  improvement  of  the  model  may  be 
made  when  satellite  infrared  data  is  available. 
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I.   INTRODUCTION 

Satellite  photographs  of  tropical  regions  have  revealed  that  much 
of  the  tropical  cloudiness  is  in  the  form  of  large,  connected  masses 
of  very  bright  clouds.   These  masses  are  composed  of  many  individual 
cumulus  clouds,  whose  extensive  cirrus  canopies  with  a  horizontal  scale 
of  several  hundred  km  or  larger  account  for  most  of  their  brightness. 
These  large,  bright  cloud  masses,  or  "cloud  clusters"  as  they  are  now 
called,  have  received  considerable  attention  in  recent  years.   Chang 
(1970)  and  Wallace  (1970,  1971),  using  time-longitude  series  of 
satellite  photographs,  found  that  many  cloud  clusters  over  the  tropical 
Pacific  are  well  organized  and  have  characteristics  resembling  those 
of  synoptic-scale  westward  propagating  wave  disturbances.    Further 
evidence  by  Reed  and  Recker  (1971)  showed  that  these  propagating  cloud 
clusters  in  the  western  North  Pacific  are  linked  to  synoptic-scale 
waves  in  the  wind  fields.   Williams  and  Gray  (1973)  and  Yanai  et  al 
(1973)  have  also  deduced  the  mean  properties  of  these  western  North 
Pacific  cloud  clusters  and  found  several  vertical  characteristics  of 
the  active  clusters  similar  to  those  found  by  Reed  and  Recker  (1971) 
and  Nitta  (1972)  for  the  wave  troughs. 

Due  to  the  poor  radiosonde  coverage  over  the  vast  tropical  oceans, 
it  is  natural  to  try  to  utilize  the  available  satellite  data  as  an 
aid  for  tropical  analysis.   Digitized  cloud  brightness  data,  acquired 
by  the  Advanced  Vidicon  Camera  System  of  the  operational  satellite 
system  have  been  available  since  1967.   Because  of  the  relationship 


between  the  bright  cloud  clusters  and  synoptic-scale  tropical  motions, 
this  data  becomes  a  potential  source  for  such  a  purpose.   In  fact, 
evidences  of  relations  between  the  satellite  brightness  data  and  other 
meteorological  parameters  have  been  found  by  several  investigators. 
The  correlation  between  brightness  and  precipitation  has  led  to  many 
attempts  to  estimate  rainfall  using  satellite  brightness  [for  example, 
see  review  by  Martin  and  Scherer  (1973)].   Wallace  (1971)  also  found 
indication  of  good  correspondence  between  area-averaged  cloud  bright- 
ness and  synoptic-scale  vertical  motion  associated  with  the  4-  to  5- 
day  tropical  disturbances  in  the  Kwajalein-Eniwetok-Ponape  triangle 
of  the  Marshall  Islands  and  the  Guam-Truk-Yap  triangle  further  west. 
He  postulated  that  it  may  be  possible  to  obtain  an  estimate  of  the  ver- 
tical motion  field  over  the  tropics  from  the  cloud  brightness  data 
alone.   If  this  turns  out  to  be  the  case,  the  satellite  cloud  brightness 
data  will  certainly  be  one  of  the  best  available  indicators  of  disturbed 
weather  in  the  tropics. 

There  are  physical  reasons  to  think  that  brightness  data  is  closely 
related  to  the  large-scale  vertical  motion  field.   The  bright  cirrus 
canopies  associated  with  the  cloud  clusters  "are  produced  by  outflow 
from  and  remnants  of  cumulonimbi.   In  general,  developing  and  conserva- 
tive clusters  maintain  their  cumulonimbi  from  a  steady  low-level  mass 
convergence;  clusters  gradually  die  when  their  low-level  mass  conver- 
gence is  eliminated"  (Williams  and  Gray,  1973).   Because  of  the  readily 
available  low-level  moisture,  the  large-scale  upward  motion  provides  a 
favorable  environment  for  the  development  and  existence  of  the  clusters. 
Furthermore,  brighter  clouds  are  usually  thicker  and  therefore  may 


indicate  larger  amount  of  latent  heat  release  (Gruber,  1974).   Because 
the  magnitude  of  the  latent  heating  is  nearly  one  order  greater  than 
the  temperature  fluctuations  at  most  levels  for  large-scale,  convectively- 
active  tropical  motions,  thermal  energy  balance  requires  that  the  large- 
scale  vertical  motion  must  be  approximately  proportional  to  heating 
(Wallace,  1971;  Holton,  1972).   Thus  one  could  expect  a  positive  corre- 
lation between  cloud  brightness  and  large-scale  vertical  motion  by 
reason  of  heating  alone. 

The  correspondence  between  brightness  and  vertical  motion  may  also 
be  extended  to  large-scale  horizontal  divergence  at  levels  where 
divergence  is  large,  since  many  observational  studies  [Wallace  (1971), 
Reed  and  Recker  (1971),  Nitta  (1972),  Williams  and  Gray  (1973),  Yanai 
et  al  (1973),  Reed  and  Johnson  (1974)]  have  produced  very  similar  pro- 
files of  the  vertical  motion  (or  divergence)  over  tropical  areas  where 
covection  is  active.   The  same  conclusion  can  also  be  inferred  from 
the  various  vertical  heating  profiles  shown  by  Hayashi  (1974)  from 
analyses  of  tropical  output  from  the  general  circulation  model  of  the 
NOAA  Geophysical  Fluid  Dynamics  Laboratory  at  Princeton.   In  addition, 
Williams  and  Gray  (1973)  have  shown  a  clear  in-phase  relation  between 
cloud  clusters  and  upper-level  divergence  in  their  composite  study. 

So  far,  the  use  of  satellite  brightness  data  as  an  aid  to  tropical 
analysis  has  been  limited  to  subjective  methods  only;  i.e.,  an  analyst 
relying  on  visual  inspection  of  the  satellite  photographs  to  modify 
conventional  weather  maps.   If  the  digitized  brightness  data  can  be  used 
to  estimate  the  large-scale  vertical  motion,  as  postulated  by  Wallace, 
and  upper-level  divergence,  as  reasoned  above,  then  there  may  be  ways 
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to  use  this  digitized  data  in  an  objective  manner.   The  purpose  of  this 
study  is  to  propose  such  a  technique  in  which  the  digitized  brightness 
will  be  used  to  objectively  diagnose  the  large-scale  flow  pattern  in 
the  tropical  upper  troposphere  over  regions  where  wind  reports  are  sparse 
or  void. 

In  Section  II,  the  design  of  the  objective  technique,  which  is 
essentially  a  barotropic  diagnostic  model  based  on  several  observed 
properties  of  the  large-scale  tropical  motions,  will  be  discussed.   In 
Section  III  the  data  used  in  the  experiments  of  the  technique  will  be 
described.   The  results  of  these  experiments  will  be  presented  in 
Section  IV. 
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II.   DESIGN  OF  THE  BASIC  MODEL 

The  objective  technique  of  using  digitized  satellite  cloud  bright- 
ness to  diagnose  large-scale  tropical  flow  is  based  on  two  primary 
suppositions: 

1)  Brightness  is  positively  correlated  with  active  convection,  and 
therefore  with  large-scale  condensation  heating,  vertical  motion  and 
upper-level  divergence. 

2)  Large-scale  tropical  flow  is  usually  quasi-barotropic  and  quasi- 
non-divergent  . 

The  model  used  is  essentially  a  linearized  barotropic  vorticity 
equation  which  may  be  applied  at  any  level  where  divergence  is  relatively 
important  in  determining  the  rotational  part  of  the  flow,  such  as  the 
maximum-divergence  level  or  the  "outflow"  level  of  the  cloud  clusters 
which  is  usually  near  200-mb  in  the  tropics.   It  is  similar  to  the  simple 
diagnostic  model  by  Holton  and  Colton  (1972)  for  the  200-mb  seasonal  mean 
circulations  during  the  northern  summer,  in  which  the  mean  motion  field 
was  successfully  reproduced  by  the  steady-state  solution  of  a  linearized 
barotropic  vorticity  equation  with  observed  mean  divergence  and  mean 
zonal  wind.   In  our  model  the  vorticity  equation  is  applied  to  a  data- 
void  region  on  a  daily  basis,  with  horizontal  divergence  replaced  by 
an  "equivalent  divergence"  field  derived  from  the  brightness  data,  basic 
zonal  wind  given  by  climatology,  and  boundary  conditions  determined  by 
the  data-rich  surrounding  areas.   Similar  to  Holton  and  Colton' s , time- 
mean  model,  a  damping  coefficient  is  introduced  to  parameterize  the 
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strong  damping  process  in  the  tropical  upper  troposphere  as  reported  hy 
several  vorticity  budget  studies  (Reed  and  Recker,  1971;  Williams  and 
Gray,  1973;  and  others).  We  hope  this  model  can  reasonably  describe  the 
rotational  part  of  the  large-scale  motion  near  the  outflow  level,  and  it 
may  be  useful  for  diagnostic  purposes  if  we  know  what  kind  of  solution 
can  best  approximate  the  actual  motion  field.   If  12-hr  or  24-hr  history 
is  available  (presumably  by  use  of  this  technique  at  an  earlier  time 
or  by  some  other  means),  extrapolation  in  time  may  be  used.   However, 
the  previous  vorticity  budget  studies  have  consistently  shown  that  the 
local  change  of  vorticity  near  the  maximum-divergence  level  is  rather 
small  when  compared  to  other  terms  in  the  vorticity  equation,  thus  a 
diagnosis  of  the  flow  at  this  level  over  the  data-void  region  may  be 
obtained  by  way  of  a  steady-state  solution  of  the  vorticity  equation. 
We  hope  that  this  solution  may,  to  some  degree,  resemble  the  actual  flow 
pattern  and  can  be  useful  for  large-scale  tropical  analysis  as  a  first 
estimate  of  the  upper-level  flow  over  a  data-void  region. 

The  design  of  the  technique  begins  from  the  linearized  barotropic 
vorticity  equation: 

^V^  +  34v^-|ii5+B|f=-(f--|)B-BV^-DV^,     (1) 

dy 


where 


The  design  of  the  technique  is  roughly  similar  to  that  described  by 
Edwards  (1973),  although  the  development  is  somewhat  incomplete  and 
there  are  also  a  few  errors  contained  in  that  paper. 
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ip  perturbation  stream  function 

C  =  V2i|j  perturbation  relative  vorticity 

u  mean  zonal  wind 

f  Coriolis  parameter 

"-■jr 

8y 

B   perturbation  divergence 

B   zonal  mean  divergence 

D   a  damping  coefficient 

x,y,t   zonal,  meridional  directions  and 
time,  respectively 

In  Equation  (1)  the  advections  due  to  mean  meridional  velocity,  mean 
vertical  velocity,  and  the  divergent  part  of  the  perturbation  meridional 
velocity  are  all  neglected  from  scale  considerations.   The  term  -  DV2iJj 
is  included  to  represent  the  strong  damping  process  in  the  upper  tropo- 
sphere as  found  by  the  budget  studies.   It  is  also  consistent  with  the 
in-phase  relationship  between  divergence  and  vorticity  near  the  outflow 
level  observed  by  Wallace  (1971)  because  such  a  relationship  requires  the 
magnitude  of  the  damping  to  be  comparable  to  the  divergence  forcing 
(Holton  and  Colton,  1972) .   This  damping  process  has  not  been  fully  under- 
stood.  It  may  be  due  to  strong  cumulus  transport  of  low-level  vorticity 
as. suggested  by  Reed  and  Recker  (1971)  and  Holton  and  Colton  (1972),  or 
it  may  also  be  due  to  some  non-linear  process.   In  their  study  Holton  and 
Colton  found  that  a  very  strong  viscous  damping  coefficient,  D  =  1.5  x 
10   sec   ,  is  necessary  in  order  to  compute  a  realistic  seasonal  mean 
stream  function  field  from  the  observed  mean  divergence.   Other  budget 
studies  (for  example,  see 
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Williams  and  Gray)  also  indicate  that  D  in  cloud-cluster  areas  is  con- 
siderably larger  than  that  may  be  expected  for  the  usual  large-scale 
damping  time  alone. 

The  corresponding  finite-difference  equation  of  (1) ,  expressed  by 
central  differencing  in  both  time  and  space (except  that  an  implicit 
scheme  is  used  for  the  damping  term  due  to  stability  requirement) , is 
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where 

T     index  for  time  step 

i,j    grid  point  indices  in  x  and  y,  respectively 

m     map  factor  (equals  to  unity  at  22. 5N  throughout 
this  study) 

d     grid  interval,  same  for  both  x  and  y  (equals 
to  2.5°  at  22. 5N  throughout  the  study) 

Equation  (2)  may  now  be  integrated  in  time  over  a  rectangular  region 
provided  that   u  ,  divergence  (B  and  B) ,  initial  conditions  and  proper 
boundary  conditions  are  given.   Two  neighboring  boundary  conditions  of 
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ty     are  required  at  both  eastern  and  western  boundaries  because  of  the 
zonal  advection  term. 

Since  we  will  be  looking  for  steady-state  solution  under  forcing, 
initial  condition  will  not  be  crucial  and  is  always  set  to  zero  for 
the  interior  in  this  study.   The  specification  of  u  posts  a  problem 
because  it  is  assumed  that  we  are  dealing  with  a  data-void  region. 
However,  it  may  be  adequate  to  use  values  known  from  monthly  or  seasonal 
mean  climatology.   The  divergence  field,  which  is  required  as  a  forcing 
function  for  (2), will  be  given  by  an  equivalent  divergence  field 
empirically  constructed  from  digitized  satellite  brightness  data.   In 
the  experiments  of  this  model  discussed  in  Section  IV,  three  forcing 
fields  for  each  day  will  be  tested  and  results  compared  with  the  avail- 
able NMC  flow  field:  a  divergence  field  kinematically  computed  from 
NMC-analyzed  winds,  an  equivalent  divergence  field  from  brightness  data, 
and  a  case  of  "zero  forcing"  where  B  =  B  =  0  . 

The  specification  of  boundary  conditions  for  ij;  will  depend  upon 
the  data  available  outside  of  the  data-void  region.   We  assume  that, 
in  most  cases,  the  data-void  region  may  be  selected  as  a  rectangular 
area  or  "box"  surrounded  on  all  four  sides  by  boxes  of  arbitrary  size 
within  which  wind  information  are  adequate  (see  Figure  1) .   These  data- 
rich  regions  may  be  areas  of  island  groups  or  edges  of  continents. 
If  stream  functions  for  each  of  these  data-rich  regions  can  be  computed 
from  the  wind  information,  the  boundary  conditions  of  the  data-void 
region  can  be  determined  because  each  of  the  four  adjacent  data-rich 
boxes  share  one  common  border  with  the  data-void  box.   There  are  many 
methods  of  computing  the  stream  function  from  winds  for  a  limited 
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region.   Among  them  the  various  versions  of  Sangster's  method,  as  dis- 
cussed by  Hawkins  and  Rosenthal  (1965),  are  all  quite  adequate  for  our 
purpose,  especially  Version  III,  which  was  found  to  give  a  stream 
function  field  closest  to  the  NMC-analyzed  field  used  in  our  test.   We 
have  also  tried  an  iterative  version  of  Sangster's  method  somewhat 
similar  to  that  described  by  Shukla  and  Saha  (1974)  but  the  difference 
between  the  iterative  and  non-iterative  version  is  usually  small,  and 
the  additional  computer  time  required  for  iteration  is  certainly  not 
justified. 

Although  theoretically  a  method  such  as  Sangster's  must  be  used 
in  the  actual  application  of  our  technique,  for  simplicity  the  boundary 
conditions  used  in  the  numerous  experiments  carried  out  in  this  study 
are  extracted  directly  from  the  available  NMC-analyzed  stream  function 
field. 
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III.   DATA  SOURCES  AND  PERIOD  OF  STUDY 

In  this  study  we  used  daily  satellite  brightness  data  for  July- 
August  1971  to  construct  the  forcing  function  for  Equation  (2) .   The 
results  are  compared  to  the  tropical  grid  wind  data  (including  stream 
functions)  analyzed  by  the  National  Meteorological  Center  (NMC) .   The 
area  selected  for  this  study  is  the  tropical  western  North  Pacific, 
from  the  equator  to  25N  and  125E  to  175E.   This  region  was  selected 
because  of  its  11  rawinsonde  stations  including  the  Caroline  and  Marshall 
Islands,  which  makes  the  NMC-analyzed  wind  data  relatively  more  reliable 
compared  to  other  areas  of  the  tropics.   Figure  2  shows  the  region  and 
the  rawinsonde  network. 

The  digitized  cloud  brightness  and  upper-level  wind  data  tapes 
were  provided  by  the  National  Center  for  Atmospheric  Research  (NCAR) . 
The  cloud  brightness  data ,  acquired  by  NCAR  from  the  National  Environ- 
mental Satellite  Service  (NESS), are  digitized  values  ranging  from  zero 
to  ten  averaged  for  5°  x  5°  latitude-longitude  squares.   The  daily 
observation  time  is  approximately  3  PM  local  time  (~  0000  GMT  in  the 
Western  Pacific) .   The  upper-level  component  winds  as  well  as  stream 
function  fields  at  the  200-mb  level  were  acquired  by  NCAR  from  the  NMC 
tropical  grid  analysis.   These  200-mb  stream  function  fields  will  be 
referred  to  as  the  "observed  fields"  when  discussing  the  experiments 
of  Section  IV.   The  wind  and  stream  function  values  are  given  at  grid 
points  approximately  coinciding  with  those  for  the  brightness  data, 
and  are  the  results  of  analyses  based  on  all  available  rawinsondes, 
aircraft  observations  and  a  few  winds  deduced  from  cloud  drift  as  seen 


from  satellite  photographs.   In  all  cases,  the  0000  GMT  wind  data  are 
considered  as  being  closest  to  the  time  of  satellite  observations.   All 
data  are  interpolated  to  give  values  at  the  2.5°  x  2.5°  latitude-longi- 
tude grid  points  used  in  this  study. 

Values  for  divergence  were  computed  kinematically  from  the  analyzed 
wind  and  found  to  be  maximum  at  200  mb.   This  is  consistent  with  observa- 
tional studies  conducted  by  Reed  and  Recker  (1971),  Wallace  (1971), 
Nitta  (1972),  Williams  and  Gray  (1973),  and  Yanai  et  al  (1973),  who  all 
found  that  a  non-divergent  level  is  usually  near  300  mb  to  400  mb  in 
the  tropics  and  the  200-mb  level  is  the  outflow  level  with  maximum 
divergence.   Therefore,  only  the  200-mb  level  is  tested  in  the  experi- 
ments of  Section  IV. 

Before  the  observed  satellite  brightness  data  acquired  from  NCAR 
can  be  used  as  equivalent  divergence  forcing  in  (2),  it  must  be  related 
to  the  200-mb  divergence  field.   This  is  accomplished  by  a  simple  scheme. 
A  proportionality  factor  K  is  computed  by  dividing  the  monthly 
standard  deviation  of  divergence,  O,.       ,  by  the  monthly  standard  devia- 
tion of  brightness,  O        .   The  equivalent  divergence  at  each  grid  point 

LB 

is  then  computed  on  a  daily  basis  by  the  following, 


(B  +  B)E  =  K  •  CB'  (3) 


where 


(B  +  B)_  =  equivalent  divergence 

b 

CB'        cloud  brightness  with  seasonal  mean 
brightness  removed 

a, . 

d  iv 
K  =  monthly  proportionality  constant 

°CB 
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When  this  equivalent  divergence  is  used  as  forcing  of  (2),  it  will  be 
referred  to  as  "brightness  forcing"  in  the  discussion  of  the  experiments 
of  Section  IV. 

The  divergence  and  brightness  used  in  this  study  have  been  analyzed 
in  a  correlation  study  by  Varona  (1974).   Weak  positive,  correlations  be- 
tween the  two  fields  were  found  in  space,  time  and  also  for  dominant 
time  scales  associated  with  synoptic  disturbances  in  the  period  of  study, 
and  the  correlation  coefficients  are  much  lower  than  those  inferred  from 
Wallace's  (1971)  wave  study.   This  may  be  due  to  the  quality  of  the  data 
used,  especially  the  divergence  which  is  kinematically  computed  from 
grid  point  winds  analyzed  by  NMC.   We  therefore  expect  that  the  results 
of  experiments  using  divergence  forcing  and  brightness  forcing  will 
generally  be  different.   It  turns  out  that  in  a  number  of  cases  the  bright- 
ness forcing  actually  produces  flow  fields  more  closely  resembling  the 
observed  field  than  those  produced  by  our  divergence  data. 
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IV.   RESULTS 

A.   TEST  OF  THE  MODEL 

As  mentioned  in  Section  II,  our  model  is  similar  to  that  used  by 
Holton  and  Colton  (1972)  for  the  200-mb  seasonal  mean  vorticity  budget 
study  during  the  northern  summer.   While  Holton  and  Colton' s  model  used 
Gaussian  elimination  after  Fourier  decomposition  in  the  zonal  direction, 
our  model  uses  a  simple  initial  value  technique  to  find  the  steady- 
state  solution.   Before  conducting  experiments  using  brightness  forcing, 
Holton  and  Colton' s  problem  is  redone  as  a  test  of  our  model.   Observed 
seasonal-mean  horizontal  divergence  and  zonal  wind  fields  for  June- 
August  1967  obtained  from  Krishnamurti' s  (1970)  data  are  used  in  Equa- 
tion (2).   Observed  stream  functions  are  specified  at  the  equator,  45N, 
87. 5W,  and  90W  as  boundary  conditions  for  this  northern  tropical  global 
band.   A  strong  damping  coefficient,  D  =  1.5  x  10   sec   ,  is  first 
used.   Holton  and  Colton  found  that  such  a  large  damping  coefficient  is 
necessary  for  the  linearized  barotropic  vorticity  equation  to  produce 
a  flow  field  which  closely  resembles  that  observed.   Figure  3a  shows 
the  observed  stream  function  field  in  this  band  as  deduced  by  Krishna- 
murti.  Figure  3b  shows  the  stream  function  field  after  72  hours  of 
integration,  which  is  considered  the  quasi-steady-state  solution  because 
further  time  integration  produces  only  slight  changes.   It  is  evident 
that  our  model,  just  as  that  of  Holton  and  Colton  (1972),  is  capable  of 
reproducing  the  major  circulation  features  forced  by  the  seasonal  mean 
divergence. 
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Holton  and  Colton  hypothesized  that  the  strong  damping  coefficient 
is  necessary  because  of  cumulus-scale  vorticity  transports  in  the 
vertical  associated  with  the  monsoonal  circulation.   However,  cumulus 
activity  is  strong  only  over  areas  with  large-scale  divergence  at  the 
200-mb  level.   Thus,  the  use  of  a  uniformly  strong  damping  coefficient 
over  the  entire  domain  seems  inconsistent  with  the  cumulus  transport 
hypothesis.   To  examine  this  premise,  a  second  test  run  of  the  200-mb 
seasonal  mean  circulation  is  carried  out  using  two  damping  coefficients. 
The  strong  damping  coefficient,  D  =  1.5  x  10   sec   ,  is  used  only  in 
the  divergence  areas  while  a  weaker  damping  coefficient,  D  =  5  x  10 
sec   ,  which  is  considered  to  be  more  realistic  for  the  large-scale  flow 
is  used  over  convergence  areas.   The  resultant  72-hr  solution  is  shown 
in  Figure  3c  and  it  is  evident  that  the  mid-oceanic  troughs  are  too 
strong.   Thus  the  necessity  of  including  a  strong  damping  coefficient 
throughout  the  domain,  in  effect,  suggests  that  some  other  viscous  pro- 
cess must  be  present  in  addition  to  the  possible  cumulus  damping. 

B.   EXPERIMENTS  WITH  BRIGHTNESS  FORCING 

In  the  daily  diagnostic  experiments  three  types  of  forcing  functions 
are  applied:   the  "divergence  forcing"  (B  +  B  =  the  kinematically  computed 
divergence),  the  "brightness  forcing"  (B  +  B  =  the  equivalent  divergence), 
and  the  "zero  forcing"  (B  =  B  =  0) .   Quasi-steady-state  solutions  after 
72  hours  of  integration  are  obtained  for  each  day  of  July  and  August  1971, 
but  only  the  results  for  an  eight  day  representative  period  are  shown 
here  (Figures  4-11).   These  are  the  first  eight  days  of  August  1971  that 
have  data  available,  namely  2,  3,  4,  5,  6,  7,  9,  and  10  August.   For 
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each  day  of  this  period  eight  computer-plotted  diagrams  of  the  perturba- 
tion stream  function  field  are  shown  in  the  following  order: 

a.  the  NMC-analyzed  field; 

b.  the  solution  for  zero-forcing  and  D  =  1.5  x  10   sec   ; 

c.  the  solution  for  divergence  forcing  and  D  =  1.5  x  10   sec   ; 

d.  the  solution  for  brightness  forcing  and  D  =  1.5  x  10   sec   ; 

e.  the  solution  for  divergence  forcing  and  D  =  5  x  10   sec   ; 

— ft  —l 

f.  the  solution  for  brightness  forcing  and  D  =  5  x  10   sec   ; 

g.  the  solution  for  divergence  forcing,  D  =  1.5  x  10   sec   ; 
and  zero  boundary  conditions  (ip  =  0  at  all  boundaries)  ;  and 

h.   the  solution  for  brightness  forcing,  D  =  1.5  x  10   sec   ,  and 
zero  boundary  conditions. 

The  stream  functions  for  2  August  are  shown  in  Figure  4.   The  solu- 
tion obtained  from  brightness  forcing  (d)  is  quite  similar  to  the  NMC- 
analyzed  field  (a),  with  some  differences  in  the  strength  of  the  major 
features  but  generally  good  agreement  between  gradients  and  flow  patterns. 
On  the  other  hand,  the  results  of  divergence  forcing  (c)  shows  a  good 
correspondence  with  the  NMC-field  over  the  eastern  half  of  the  test 
region  but  it  contains  a  ficticious  trough  in  the  northwestern  "high" 
area.   The  zero-forcing  solution  (b)  least  resembles  the  NMC  field, 
especially  in  the  north-central  sector  where  the  flow  displays  excessive 
zonal  orientation.   As  expected  the  flow  fields  are  too  strong  for  both 
the  divergence  and  brightness  forcing  when  the  damping  coefficient  is 
reduced  to  D  =  5  x  10   sec    (e  and  f ) .   The  last  two  diagrams  illustrate 
the  results  of  the  zero  boundary  condition  cases.   Here  the  brightness 
forcing  (h)  produces  a  pattern  that  better  defines  the  distribution  of 
positive  and  negative  vorticity  areas  than  does  the  divergence  forcing 
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(g)  which  is  consistent  with  the  results  obtained  when  realistic 
boundary  conditions  are  given. 

Figure  5  shows  the  stream  functions  for  3  August.   It  can  be  seen 
that  brightness  forcing  produces  a  flow  field  in  excellent  agreement 
with  the  NMC  field  for  the  entire  region.   The  field  produced  by 
divergence  suffers  mainly  from  a  fictious  trough  which  is  present  in 
the  western  "high"  area.   The  zero-forcing  solution  again  least 
resembles  the  NMC-field.   In  the  eastern  central  portion  the  flow  is 
both  weak  and  more  zonal,  and  the  cyclonic  vorticity  also  penetrates 
into  the  western  area  although  not  as  strongly  as  the  trough  shown  in 
the  divergence-forced  field.  The  results  for  smaller  damping  are  again 
too  strong.   This  is  consistently  true  for  all  days  tested  and  we  will 
not  make  further  special  reference  to  it.   The  value  D  =  1.5  x  10  sec 
is  thus  considered  an  appropriate  damping  coefficient  to  be  used  in 
our  model.   The  results  for  the  zero  boundary  conditions,  as  expected, 
show  that  the  brightness  forcing  better  delineates  the  regions  of  nega- 
tive and  positive  vorticity  than  does  the  divergence  forcing. 

The  results  for  4  August,  as  shown  in  Figure  6,  differ  from  those 
obtained  on  the  two  previous  days.   The  field  produced  by  brightness 
forcing  does  not  agree  as  closely  with  the  NMC  field  as  it  did  on  2  and 
3  August.  Both  the  western  and  eastern  "cells"  are  somewhat  under- 
developed.  However,  the  orientation  of  the  eastern  trough  is  better 
depicted  by  this  field  than  by  the  field  produced  by  the  divergence  or 
zero  forcing.   On  the  other  hand,  the  western  anticyclonic  circulation 
pattern  produced  by  the  divergence  forcing  more  closely  resembles  that 
of  the  NMC  field.   The  accurate  orientation  of  the  eastern  trough 
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produced  by  the  brightness  forcing  is  clearly  depicted  in  the  results 
for  zero  boundary  conditions.   This  is  seen  in  Figure  6h,  which  clearly 
defines  a  zone  of  negative  stream  function  values  sloping  from  the 
northeast  to  the  southwest.   This  orientation  does  not  appear  in  the 
field  produced  by  the  divergence  forcing  with  zero  boundary  conditions. 

Results  for  5  August  are  shown  in  Figure  7.   The  major  circulation 
features  produced  by  the  brightness  forcing  do  not  develop  meridionally 
into  the  south  central  portion  as  those  of  the  NMC  field.   This  parti- 
cularly true  of  the  southern  extension  of  the  -1  contour  line  of  the 
eastern  trough.   Somewhat  stronger  development  of  this  feature  is  pro- 
duced by  the  divergence  forcing.   This  feature  is  least  developed  in 
the  solution  obtained  by  zero-forcing,  thus  on  this  day  the  brightness 
forcing  result  may  still  be  considered  as  better  than  the  zero-forcing 
case.   The  zero  boundary  condition  results  show  that  although  the 
divergence  forcing  case  defines  the  orientation  of  the  trough  to  a 
certain  degree  of  accuracy,  the  brightness  forcing  case  depicts  the 
position  and  strength  of  the  trough  much  better. 

Figure  8  shows  that  for  6  August  brightness  forcing  produces  areas 
of  anticyclonic  vorticity  over  the  western  half  and  southeastern  corner 
of  the  test  region,  both  of  which  grossly  resemble  the  patterns  shown 
by  the  NMC-field.   However,  the  northeastern  cyclonic  circulation  is 
overdeveloped  by  the  brightness  forcing.   It  can  be  seen  that  this 
northeastern  quarter  of  the  test  area  shown  in  the  zero-forcing  solution 
has  closer  agreement  with  the  NMC-field  than  does  the  divergence  or 
the  brightness  forcings,  but  the  overall  resemblance  is  still  best  for 
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the  field  produced  by  brightness  forcing  because  of  the  two  anticyclonic 
areas.   When  zero  boundary  conditions  are  used,  the  brightness  forcing 
solution  is  again  considerably  more  accurate  in  specifying  the  negative 
and  positive  vorticity  regions  than  the  results  of  divergence  forcing, 
although  it  is  not  as  representative  of  the  major  features  of  the  NMC- 
field  as  is  the  case  for  several  previous  days. 

The  results  for  7  August  are  shown  in  Figure  9.   The  correspondence 
between  the  solution  obtained  with  brightness  forcing  and  the  NMC-field 
appear  to  be  quite  poor.   However,  closer  examination  reveals  that  the 
slope  of  the  -1  contour  produced  by  brightness  forcing  closely  approxi- 
mates that  of  the  zero  contour  of  the  NMC-field,  which  is  orientated 
NW-SE  near  the  central  western  portion.   Therefore,  even  though  the 
appearance  of  the  two  fields  differs  considerably,  the  brightness  forc- 
ing solution  still  gives  some  useful  information  on  the  NW-SE  flow 
direction  in  the  central  part  of  the  test  region,  which  is  not  the  case 
for  the  fields  obtained  from  divergence  or  zero-forcing.   This  infor- 
mation on  flow  direction  is  also  evident  in  the  -1  contour  line  of  the 
brightness  forcing  solution  with  zero  boundary  conditions. 

Figure  10  shows  the  stream  functions  for  9  August.   On  this  day  all 
results  are  considered  poor  and  is  clearly  a  failing  case  of  our  model. 

The  results  for  10  August  as  shown  in  Figure  11  illustrates  another 
unsuccessful  case.   The  solution  obtained  from  brightness  forcing, 
however,  still  resembles  the  NMC-field  better  than  the  other  forcing 
solutions.   The  patterns  produced  by  zero  boundary  conditions  are 
obviously  poor  in  terms  of  representing  the  major  circulation  features 
of  the  NMC-field. 
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In  summary,  the  above  results  show  that  the  solutions  obtained  with 
the  strong  damping  coefficient,  D  =  1.5  x  10   sec   ,  are  generally 
more  realistic  than  those  obtained  with  D  =  5  x  10  .  sec   .   In  the  8 
days  discussed,  brightness  forcing  produced  stream  functions  quite 
similar  to  the  NMC-field  on  two  days  (2  and  3  August)  and  gives  useful 
information  on  the  flow  patterns  on  four  days  (4,  5,  6  and  10  August). 
It  is  also  interesting  to  note  that  on  the  failing  day  of  August  9, 
both  the  computed  divergence  and  the  brightness  values  are  generally 
small.   On  this  day  few  values  of  brightness  data  in  the  test  area 
exceed  7,  which  is  the"lowest  brightness"  day  of  the  8-day  period. 

In  addition  to  the  solutions  presented  here,  results  were  obtained 
for  all  days  in  the  months  of  July  and  August  for  which  data  are 
available.   Overall,  for  approximately  25%  of  the  tested  days,  brightness 
forcing  produces  quite  realistic  stream  functions  resembling  the  NMC- 
field,  and  in  another  30%  of  the  days  it  gives  some  useful  information 
about  the  flow  field.   These  results  suggest  that  if  no  wind  informa- 
tion is  available  over  a  tropical  region,  a  first  guess  of  the  200-mb 
flow  field  may  be  obtained  from  digitized  satellite  brightness  data 
using  a  simple  model  such  as  that  presented  here. 

Since  the  zero  forcing  results  occasionally  also  give  a  rough  sketch 
of  the  flow  pattern,  even  steady-state  solutions  of  the  barotropic 
vorticity  equation  from  boundary  information  alone  may  sometimes  provide 
a  guess  of  the  flow  which  is  better  than  nothing.   This  is  especially 
true  when  divergence  is  weak  as  indicated  by  low  brightness  values. 
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C.   EXPERIMENTS  WITH  FORCING  MODIFIED  BY  TIME  INTERPOLATED  INFRARED  DATA 

Since  the  satellite  digitized  brightness  does  not  supply  direct  in- 
formation on  the  heights  of  the  cloud  clusters  observed,  the  possible 
existence  of  low-level  high  brightness  values  may  distort  the  upper- 
level  flow  pattern  obtained  from  the  solution  of  Equation  (2)  under 
brightness  forcing.   If  satellite  infrared  data  (IR)  were  available  for 
the  same  time  and  grid  points  as  satellite  cloud  brightness  data,  the 
values  of  brightness  could  be  modified  according  to  the  height  of  the 
cloud  tops,  which  can  be  inferred  from  the  IR  data.   This  is  because 
the  condensation  heating  at  upper  levels  is  obviously  related  to  the 
vertical  extent  of  cumulus  development.   If  accurate  cloud-top  tempera- 
tures and  cloud-top  heights  can  be  determined,  there  may  be  ways  to 
incorporate  this  information  into  the  brightness  data  in  order  to  dis- 
tinguish the  effects  due  to  various  vertical  developments  of  cumulus 
convection.   For  example,  since  higher  cloud  tops  are  colder,  we  may 
assume  that  by  giving  higher  weights  to  brightness  values  coincident 
with  low  IR  readings  the  resultant  field  obtained  by  brightness  forcing 
of  Equation  (2)  will  better  approximate  the  actual  field.   In  this 
way  the  high  brightness  due  to  solar  reflections  at  the  surface  would 
also  be  completely  eliminated  as  far  as  the  contribution  to  the 
divergence  at  the  upper  level  is  concerned. 

In  1971,  area-averaged,  once-daily  IR  data  from  the  NOAA-1  satellite 
over  the  5°  x  5°  latitude-longitude  grids  of  the  test  region  are  avail- 
able from  NESS  for  the  period  19  April  -  10  July.   Unfortunately,  the 
observation  time  is  approximately  3  am  local  time  each  day  (or  ~  1200 
GMT  in  the  western  Pacific),  which  is  12  hours  different  from  that  of 
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the  brightness  data.   However,  although  the  life  cycle  of  an  individual 
cumulus  element  may  be  ~  0  (hr),  an  organized  cloud  cluster  with  a  scale 
on  the  order  of  the  5°  x  5°  grid  squares  usually  lasts  sveral  days  or 
longer.   Thus  we  decided  to  use  the  simple  average  of  two  consecutive 
IR  observations  to  approximate  the  IR  information  at  the  brightness 
observation  time.   This  interpolated  data  is  then  used  to  test  the  idea 
of  enhancing  the  brightness  forcing.   Of  course,  this  interpolation  pro- 
cess may  introduce  large  amounts  of  error,  especially  since  cloud 
clusters  are  usually  in  a  transient  state  with  westward  movement  com- 
parable to  the  speed  of  the  zonal  flow  at  low  levels  (Chang,  1970). 
Nevertheless,  this  crude  representation  of  the  IR  data  may  provide  us 
with  some  clue  as  to  whether  improvement  can  be  reasonably  expected  in 
cases  when  actual  IR  data  is  available  for  use. 

The  NOAA-1  IR  data  period  overlaps  our  experimental  period  only 
for  the  first  ten  days  of  July  1971;  and,  due  to  missing  data  in  the 
IR  tapes,  only  IR  data  for  4-5  July  and  7-8  July  can  be  used  in  our 
averaging  procedure.   Thus  the  "0000Z  IR  fields"  are  constructed  for 
only  two  days:  5  July,  which  is  the  average  of  IR  data  of  4  July  and 
5  July;  and  8  July,  which  is  the  average  of  7  and  8  July.   As  done  pre- 
viously for  the  brightness  and  NMC  data,  these  averaged  IR  data  are 
further  interpolated  to  give  values  at  the  2.5°  x  2.5°  grid  points. 
A  quite  simple  weighting  scheme  is  then  used  to  modify  the  brightness 
values  by  the  averaged  IR  field.   This  scheme,  which  is  based  on  the 
facts  that  black-body  radiation  is  proportional  to  the  forth  power  of 
the  temperature  while  temperature  varies  almost  linearly  with  height, 
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is  described  by  the  following  equation: 

CBM  =  (C  -  4^IR)  CB  ,  (4) 

[CB  =0   if  right  hand  side  of  (4)  is  negative] 

where 

CB,,   IR  modified  brightness 
M 

CB    original  brightness 
IR    averaged  IR  data 
C     a  parameter 

The  parameter  C  is  meant  to  represent  the  mean  sea-surface  temperature 
which  may  have  different  values  for  different  latitudes  and  is  the 
fourth  root  of  a  high  IR  radiation  reading  at  each  latitude.   For  our 
simple  testing,   C   is  determined  at  each  latitude  in  the  following 
manner :  / 

1)  The  maximum  IR  value  before  averaging  is  determined  for  each 
day  (4,  5,  7,  and  8  July); 

2)  The  minimum  of  these  four  maximum  values  is  taken  as  the  repre- 
sentative sea  surface  IR  radiation;  and 

3)  C  is  then  set  to  equal  to  the  fourth  root  of  this  IR  radiation. 
The  brightness  data  as  modified  by  (4)  is  now  transformed  to  the 

equivalent  divergence  field  by  the  method  described  in  Section  III. 
Equation  (2),  with  D  =  1.5  x  10   sec   ,  is  then  used  to  obtain  72-hour 
quasi-steady-state  solutions  for  5  and  8  July  1971,  and  compared  with 
the  results  obtained  previously.   These  are  illustrated  in  Figures  12 
and  13  in  which  the  following  format  is  used: 
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a.  the  NMC-analyzed  field, 

b.  the  solution  for  divergence  forcing, 

c.  the  solution  for  brightness  forcing,  and 

d.  the  solution  for  IR-modified  brightness  forcing. 

Figure  12  shows  the  results  for  5  July  1971.   It  can  be  readily  seen 
that  the  gradients  and  pattern  of  the  northwestern  "high"  area  produced 
by  the  brightness  forcing  is  improved  by  the  IR  modification  and 
resembles  the  NMC-field  better.   On  8  July  1971,  which  is  shown  in 
Figure  13,  significant  improvement  can  again  be  found.   In  the  center 
and  center-east  sections,  the  gradients  produced  by  the  IR  modification 
is  evidently  closer  to  those  of  the  NMC  field  than  the  unmodified  forc- 
ing result.   Thus  our  simple  test  suggests  that  there  is  great  potential 
for  the  use  of  IR  data  to  enhance  the  digitized  brightness.   The  prospect 
in  the  future  should  be  quite  bright  when  concurrently-available  IR 
data  are  used  with  brightness  data,  and  when  a  better  modification 
method  is  designed  to  replace  (4). 
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V.   SUMMARY  AND  CONCLUSION 

In  this  study  we  proposed  a  simple  objective  technique  which  may 
be  used  to  obtain  large-scale  flow  information  over  tropical  regions 
where  wind  reports  are  sparse  or  void.   The  technique  is  basically 
a  diagnostic  model  consisting  of  a  linearized  barotropic  vorticity 
equation.   Primary  assumptions  of  the  technique  are:   1)  satellite 
brightness  data  is  positively  correlated  with  horizontal  divergence; 
2)  tropical  flow  is  quasi-barotropic  and  quasi-non-divergent.   Partial 
observational  justification  for  the  first  assumption  has  been  provided 
by  Wallace  (1971),  and  the  second  one  is  accepted  from  the  results 
of  scale  analysis  (Holton,  1972).   Another  assumption  used  in  the 
present  technique  is  that,  at  levels  where  divergence  is  strong,  the 
rotational  part  of  the  large-scale  flow  is  largely  determined  by  the 
balance  between  divergence  forcing,  strong  viscous  damping  and  advec- 
tion.   This  "steady  state"  assumption  is  based  on  the  results  of 
various  budget  studies  (for  example,  Reed  and  Recker,  1971;  and 
Williams  and  Gray,  1973). 

The  technique  employs  an  initial-value  approach  to  a  finite-differ- 
ence form  of  the  linearized  barotropic  vorticity  equation  applied  at  a 
strong-divergence  level  (usually  near  200-mb) .   The  strong  damping 
process  is  parameterized  using  a  damping  coefficient.   Three  types  of 
information  are  needed  to  obtain  the  quasi-steady-state  solution  of 
this  equation,  which  is  considered  to  be  an  estimate  of  the  actual  flow 
field.   These  are,  1)  appropriate  boundary  conditions,  2)  the  mean 
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zonal  flow,  and  3)  a  divergence  forcing  function.   It  is  suggested  that 
Sangster's  (1960)  method  for  computing  stream  functions  from  winds  can 
be  used  for  obtaining  boundary  conditions  from  the  wind  reports  in 
areas  surrounding  the  data-void  region.   The  mean  zonal  winds  are 
supposedly  available  from  monthly-mean  climatology.   The  forcing  func- 
tion is  an  equivalent  divergence  field  empirically  constructed  from  the 
digitized  satellite  brightness  data. 

The  model  is  tested  at  the  200-mb  level  for  each  day  of  July  and 
August  1971,  over  a  region  of  the  western  North  Pacific  bounded  by  the 
equator,  25N,  125E,  and  175E.   The  data  base  is  formed  by  the  daily 
digitized  cloud  brightness  data  from  NESS,  and  the  analyzed  upper-level 
winds  and  stream  functions  from  NMC.   Both  of  these  are  available  at 
5°  x  5°  latitude-longitude  grid  points  and  are  packed  into  computer 
tapes  by  NCAR.   A  daily  divergence  field  is  computed  kinematically  over 
the  test  region  using  the  analyzed  wind  field  based  on  0000  GMT  wind 
observations.   Comparisons  are  then  made  among  the  72-hour  quasi-steady- 
state  solutions  of  the  model  as  it  responds  to  forcing  by  the  kinematic 
divergence  field  (divergence  forcing) ,  an  equivalent  divergence  field 
constructed  from  brightness  data  (brightness  forcing),  and  a  case  in 
which  no  forcing  is  applied  (zero  forcing) .   Results  for  a  representa- 
tive period,  the  first  8  days  of  August  with  available  data,  are 
presented  and  discussed. 

In  approximately  25%  of  the  days  of  the  entire  two-month  period  of 
experiments,  the  solutions  produced  from  brightness  forcing  are  found 
to  accurately  reproduce  the  major  features  of  the  perturbation  stream 
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function  fields  analyzed  by  NMC.   In  another  30%  of  the  days  these  solu- 
tions provide  some  valuable  information  about  the  flow  patterns  which 
would  be  difficult  to  obtain  from  boundary  information  alone.   In  some 
cases,  especially  when  the  divergence  and  brightness  values  are  low, 
the  solutions  of  zero  forcing  also  resemble  the  NMC  field  to  a  certain 
degree.   The  results  of  these  experiments  thus  suggest  that  for  data- 
sparse  tropical  regions,  this  model  may  be  useful  in  obtaining  a  first 
estimate  of  the  actual  large-scale  flow  field.   This  estimate  may  then 
be  used  to  supplement  whatever  other  information  is  available  to  improve 
the  analysis  in  these  regions. 

The  technique  proposed  here  may  be  improved  in  the  future  when 
better  data  and  better  understanding  of  the  various  physical  processes 
of  the  large-scale  tropical  flow  are  available.   For  example,  it  is 
shown  in  two  cases  of  this  study  that  even  with  crude  representations  of 
the  satellite  IR  data,  the  brightness  forcing  can  be  enhanced  to  produce 
more  realistic  flow  patterns.   Better  understanding  of  the  physical 
processes  should  enable  us  to  improve  the  parameterization  of  the 
strong  damping  and  the  handling  of  the  local  time  change  term.   Both 
better  data  and  knowledge  of  the  physical  processes  would  also  facili- 
tate our  understanding  of  the  relationships  between  satellite  cloud 
data,  condensation  heating  and  the  resultant  divergence  and  vertical 
motion.   Therefore,  the  scheme  for  specifying  the  forcing  field  from 
the  various  satellite  data  may  also  be  greatly  improved. 
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Figure  4.   Perturbation  stream  functions  at  200-mb  for  2  August  1971 
(contour  interval  is  5  x  10  m  sec  ).   Fields  shown  are: 

a.  The  NMC-analyzed  field;  and  the  72-hr  solutions  resulting  from 

b.  zero  forcing  with  D  =  1.5  x  10   sec   , 

c.  divergence  forcing  with  D  =  1.5  x  10   sec   , 

d.  brightness  forcing  with  D  =  1.5  x  10   sec    , 

e.  divergence  forcing  with  D  =  5  x  10   sec   , 

f.  brightness  forcing  with  D  =  5  x  10   sec   , 

g.  divergence  forcing  with  zero  boundary  conditions  and 
D  =  1.5  x  10   sec   ,  and 

h.   brightness  forcing  with  zero  boundary  conditions  and 
D  =  1.5  x  10~5  sec"1  . 
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